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Abstract 

Kravchuk orthogonal functions Kravchuk polynomials multiplied by the square root 
of the weight function — simplify the inversion algorithm for the analysis of discrete, finite 
signals in harmonic oscillator components. They can be regarded as the best approximation 
set. As the number of sampling points increases, the Kravchuk expansion becomes the 
standard oscillator expansion. 


1 Introduction 

In a harmonic oscillator environment, such as Fourier optics in a multimodal parabolic index-profile 
fiber, sampling on a finite set of discrete observation points reconstructs the wavefunction 
through partial wave synthesis. For the harmonic oscillator eigenfunctions, one must invert a 
nondiagonal matrix with the dimension of the number of data. 

We show that Kravchuk orthogonal functions optimize the algorithm for the expansion coef- 
fients, because the matrix is already diagonal. 

• Kravchuk functions [1, 2] are solutions of the difference analogue of the Schrodinger equation 
describing a discrete harmonic oscillator system. 

• Kravchuk functions have a well-defined analytical structure inside the measurement interval. 

• Kravchuk functions become the standard oscillator wavefunctions, as the number of sampling 
points increases. 

i 

This contribution is a resume of Ref. [3]. 
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2 Harmonic oscillator expansions over a lattice of sam- 
pling points 


The standard harmonic oscillator eigenfuctions are 

1 


iMO = 


\J \/7r2 n n! 


H n (Oe~ e/ \ n = 0,1,2,..., 


where // n (£) are the Hermite polynomials, £ = yjmu/hx, m is oscillator mass, u is oscillation 
frequency, and the position coordinate is x. This function set is orthonormal under the £ 2 (5R) 
inner product: 

too f 1, when m = n, 

(^mi — m,n when TU ^ Ti. 

Thus, an arbitrary function /(£) € £ 2 (&) can be approximated in the norm as 

OO 

/(0 = £ 

n=0 

where the expansion coefficients {c n }^L 0 are determined by 

/ OO 

diMO /( 0 - 

-OO 

When the TV + 1 values {/(&)}f =0 of function /(£) are sampled on the points 
£o = ~\Nh, ... £j = (~\N + j) h,... 6v = \Nh , 


then 


= j =0,1,... N. 

n= 0 

The task to determine the iV + 1 coefficients {c (N) }^ =0 is formulated in matrix form as 

f = *W c ( N) , 

where 


( 2 . 1 ) 


II 

m o) i 
/«.) 

II 




/«») J 


«B° / 


and 


y(N) _ 


/ ^o(£o) </>o(£i) ••• M(n) \ 

V>i(&») <M6) ••• V’iUn) 


\ tM£o) Vw(6) ••• ^n{Cn) J 

This is the TV x JV matrix that has to be inverted to find the coefficients in (2.1). 
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3 Kravchuk functions are difference analogs of the osci- 
llator eigenfunctions 

Kravchuk polynomials k^\x, N) are: 

• polynomials of degree 0 < n < N, 

• in the variable x £ [0, V], 

• of the parameter 0 < p < 1. 

These polynomials are related to the binomial distribution of probability theory [4, 5], They form 
an orthogonal set 

EeU)&i,N)kMu,N) = 

3=0 

with respect to a discrete binomial weight function 

<>(*) = cSp*(i -p) n -. 

Kravchuk functions are defined as 

<t>n\ x , N) = d ~ 1 k (p) (Np + x , N) p 1/2 (x + Np), 

0 <n<N, —Np < x < (l — p)N 

(cf. definition of the Hermite functions in Ref. [6]). They obey the three-term recurrence relation 

[x - n - p(N - 2 n)]k { n p) (x, N) = (n + 1 N) + p(l - p)(N - n + 1)^, (x, N), 

and satisfy the equation 

with the difference Hamiltonian 

H (W) (x) = (1 - 2 p)x + 2p(l - p)N + \ - yjp{ 1 -p)[a(x) e~ a * -f a(x + 1) e a *], 

namely 

[(1 -2p)x-n + 2p{l -p)N)(t>W{x,N) = yjp(l - p) [a{x)(t>^\x - 1, N) + a(x + 1 )<$\x + 1, N)]. 

The oscillator equation of motion in the Schrodinger representation [7] is [//,[//, x]] = (hu>) 2 x. 
The difference analogue of this relation satisfied by this Hamiltonian is [1] 

[Hl' v >( I ),[Hl"l(x),x]] = x. 

Finally, the limit N — * oo of Kravchuk functions is 

lim /i _1/2 4 p) (/i _1 C^V) = V’n(0- 

fy — ► oo 
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The set of Kravchuk functions *«(*,*), n = 0,1,...,JV forms a basis for irreducible represen- 
tations of the rotation group 50(3) [8), corresponding to the eigenvalues t = jW of' the invariant 
Casimir operator; the eigenvalues of generator J z are the integer m - n 2 N - n t. 1 e 
representations corresponding to different values of the parameter p turn out to be umtanly equi- 
valent [1], so it is sufficient to consider a set of functions fi^(x,N) with some fixed value of this 
parameter. It is convenient to choose the value p = |, since these Kravchuk functions have definite 
parity with respect to reflections of x, 

We thus use henceforth the symmetric Kravchuk functions 

<f> n {x,N) = 2 n ~ N/2 k n {x+ l -N,N) 


4 Finite approximation by Kravchuk, functions 

A function /(£) that ‘lives in a harmonic oscillator environment’, of which the values on N + 1 
equidistant points are known, is meaningfully expanded in symmetric Kravchuk functions as 

Jfte) = 4? IX W) *.&/*.*). ; =0,1,... IV. 

I V fl 71 = 0 } 

To find the expansion coefficients {4 N) }n= cm we multiply the above equation by <f> m {£j/h , N) and 
sum over the sample points : 

4 ") = Mii/k, N)f(ii). 

j = 0 

We thus have only to multiply the sampled values /(&) by the (numerically calculated) values 
of the Kravchuk functions at the points x 3 = &/fe, for n = 0,1,..., A, to find the expansion 
coefficients. No matrix inversion is necessary. 

The sum N 

/({, N) = -L £ 4 N) *.«/*. ">• > = 0,1 iV, 

V "> n = 0 

interpolates the original function defined on discrete points to the interval [6, 6v] and is a finite 
approximation to the square-integrable function /(£)■ This approximantjs finite because for 
any fixed A’ it has a finite support {-h - h~\h + h" 1 ) with h = y/2/N. When A grows, 
the approximation to /({) becomes better. The time evolution of the approximating function 
multiplies each 4> n (tj,N) by the usual time dependence exp (-iE n t/h), with the equally spaced 

energy eigenvalues E n = hu(n + |). 
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5 Position and momentum functions 

The canonical vector basis of position functions 



is interpolated to 

f((i) = E ll N> JV) = /}">, i = o, 1, . . . N. 

k=0 

These basis functions can be expanded in terms of Kravchuk functions as 

A Mi,N) = -4 E C Mii/h, n , = 0, 1, . . . N, 

V ft n=0 

where the coefficients are __ 

i iff = </hM(hlh,N) t 

for jfc = 0, 1, . . . [N/2] and continuous £. These functions are the localized states of the discrete 
oscillator. 

Momentum basis functions are defined in the same way, because Kravchuk functions are 
self-reproducing under the discrete Fourier transformation [2], i.e. 

ML A0 = £* n ML/h, N) <t> n (UK n). 

n = 0 

Acknowledgments 

We are grateful to Mesuma K. Atakishiyeva and Guillermo Krotzsch for supporting computation 
and graphics. One of us (N.M.A.) would like to thank Institute de Investigaciones en Matematicas 
Aplicadas y en Sistemas, UN AM (Cuernavaca) for hospitality. This work is partially supported 
by the DGAPA Project N 104293. 

References 

[1] N.M.Atakishiyev and S.K.Suslov, Theor. Math. Phys. 85 , 1055 (1991). 

[2] R. Askey, N.M. Atakishiyev, and S.K. Suslov, Proceedings of the XV Workshop on High 
Energy Physics and Field Theory (Protvino, Russia, 6-10 July 1992), to appear. 

[3] N.M. Atakishiyev and K.B. Wolf, Rev. Mex. Fis. 40, 366 (1994). 

[4] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. Tricomi, Higher Transcendental Functions , 
Vol.2 (McGraw-Hill, New York, 1953). 


407 



[5] G. Szego, Orthogonal Polynomials (Amer. Math. Soc., New York, 1959). 

[6] N. Wiener, The Fourier Integral and Certain of Its Applications (Dover, New York, 1933). 

[7] L.D. Landau and E.M. Lifshitz, Quantum Mechanics: Non- Relativistic Theory (3 rd ed., Per- 
gamon Press, Oxford, 1963). 

[8] I.M. Gel’fand, R.A. Minlos, and Z.Ya. Shapiro, Representations of the Rotation and Lorentz 
Groups and Their Applications (Pergamon Press, Oxford, 1963). 


408 



